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Abstract 

The different magnetic phases of the intermetallic compound NdsGes are studied in terms of the 
specific heat, in a broad range of temperatures (350 mK-140 K) and magnetic fields (up to 40 kOe). 
The expected and terms are not found in the antiferromagnetic (AFM) and ferromagnetic 
(FM) phases respectively, but a gapped contribution that originates from a mixture of AFM 
and FM interactions in different dimensionalities under a large magnetocrystalline anisotropy, is 
present in both. An almost identical Schottky anomaly, that arises from the hyperfine splitting of 
the nuclear levels of the Nd^"*" ions, is observed in both phases, which leads us to state that the 
magnetic-field induced transition AFM—>-FM that the system experiments below 26 K consists in 
the flip of the magnetic moments of the Nd ions, conserving the average local moment. 
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I. INTRODUCTION 


In the last years the binary intermetallic compound NdsGes has been the object of dif¬ 
ferent studies. The interest in this system comes from the fact that many of its physical 
properties present abrupt changes when the transition between its two magnetic phases oc¬ 
curs. The first investigations of the magnetic properties showed that this material orders 
ferrimagnetically in zero applied magnetic field at the Neel temperature, T^r 50 K, and 
possesses a remanent moment at 4.2 K-. Later on, neutron diffraction experiments sug¬ 
gested that when the system is cooled in zero applied magnetic field an antiferromagnetic 
(AFM) state is established below T^r, although the hysteresis loop typical of a hard magnetic 
material was found at 4.2 K. This fact indicated a magnetic-field-induced phase transition 
to a ferromagnetic (FM) stated. In recent years, the system has regained attention fostered 
by the exploration of the temperature dependence of the irreversibility of this transition 
below Tt = 26 K-. Some papers have been published on this compound showing rich phe¬ 
nomena in magnetostriction^, electric resistance, specific heat, spontaneous magnetic phase 
transitions^, and more recently optical properties and electronic structure^. 

The intermetallic alloy NdsGea belongs to the family RsGea, with R = rare earth, that 
crystallizes in the hexagonal MusSis-type structure {PG^/mcm, space group No. 193). The 
structure contains two formula units per unit cell, in which R atoms occupy the two non¬ 
equivalent crystallographic sites Ad and 6g, while Ge atoms occupy the 6g site^. Although it 
is still unclear, the most accepted magnetic structure, which has been derived from neutron 
diffraction experiments^'^'^, consists in a complex double sheet, with the magnetic moments 
of the Nd ions located in the 6g position oriented parallel to the c-axis, and the moments of 
the Nd ions in the Ad position oriented along the c-axis with a deviation angle of 31° and a 
propagation vector k = (0.25 0 0); the z component changes sign every two successive (110) 
planes. 


f.y 


and 


The temperature variation of the specific heat in zero field was measured in Ref. 

10 to examine the magnetic phase transitions. Both works observed a hump around 50 K 
but no anomaly around 26 K, which was attributed to the existence of a spin-glass state, 
because of the occurrence of similar features in the Cp{T) curve in other well-known spin- 
glass systems. In this paper we will conduct a detailed study of the specific heat, as a 
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function of the temperature and the magnetic held, seeking a better understanding of the 
properties of the different magnetic phases. 


II. METHODS 

Polycrystalline ingots were prepared by arc-melting the constituting 99.9%-pure Nd and 
99.999%-pure Ge elements under high-purity argon atmosphere. The compounds were found 
to be single-phase by powder X-ray diffraction. Single crystals were grown by the Czochralski 
method from single-phase polycrystalline samples using a tri-arc furnace. It should be noted 
that it is difficult to grow a large single crystal of NdsGes because the grown crystal ingots 
tend to have small single-crystalline grains. The sample was cut from one ingot into a 
rectangular shape (1 x 1.5 x 2 mm^) and annealed at 300°G for 24 h in an evacuated 
quartz tube. The crystal orientation was determined by the back-rehection Lane method. 
Measurements of the specihc heat in the temperature range from 350 mK to 300 K and 
magnetic helds up to 40 kOe were made using the heat pulse-relaxation method with the 
heat capacity option of the PPMS® system, produced by Quantum Design®. 


III. RESULTS AND DISCUSSION 

Figure [U shows the specific heat, G, of our NdsGes single crystal (triangles) and a poly- 
crystaline sample of the nonmagnetic isostructural compound LasGes (diamonds; extracted 
form Joshi et alM-) as a function of temperature. The specific heat of the latter will be used 
as a blank and follows the expected monotonic behavior from the electronic and phononic 
contributions!^, while the specific heat of NdsGes presents a hump around ~ 50 K and 
thenceforth tends progressively to the Dulong-Petit limit. 

The dependence of the specific heat of NdsGes with the temperature was investigated as 
follows. Following a zero-field cooling process (ZFG) the specific heat was measured from 
140 K down to 350 mK. During this process we observed the expected hump associated 
to the paramagnetic-AFM transition at 50 K. The next step was to measure the system in 
the FM state. It is known that cooling the NdsGes compound below 26 K with an applied 
magnetic field larger than 5 kOe leaves the system in the saturated FM stated. To ensure 
this fact we measured the specific heat following a field cooling process (FG) in an applied 
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FIG. 1. (Color online) Zero-field temperature dependence of the specific heat of a single crystal of 
NdsGes (triangles) and a polycrystal of LasGea (diamonds; extracted from Joshi et a/.-ii). The lines 
joining the data points are guides to the eye. The vertical dotted line indicates the temperature 
(50 K) at which a hump is observed for NdsGes. The horizontal solid line indicates the Dulong-Petit 
limit of the specific heat. 

field of 15 kOe, large enough for our purpose. Nonetheless our goal was to compare the 
magnetic contribution of the different magnetic phases to the specihc heat, and the applied 
magnetic held could play an undesired role. Therefore, after measuring the FC process we 
set the magnetic held to zero and measured the specihc heat of the ferromagnetic remanent 
state (FMrem) as we increased the temperature. From magnetization experiments it is known 
that well below 26 K the FM^m and the FMpc are magnetically equivalent^. We observed 
that this equivalence is also present in terms of specihc heat, as it is shown in Fig. |2l In this 
hgure the three data sets are plotted together showing the AFM curve and how the FMrem 
and FMpc curves superimpose. Consequently, from now on we will rename the FMrem state 
as FM state in this temperature region. 

The specihc heat can be assumed to be made up of four independent contributions, 

C{T) = CeKT) + CUT) + ayp(r) + C'rnag(T). (1) 

The contribution from phonons, Ciat, can be subtracted using the specihc heat of the 
nonmagnetic isostructural compound LasGes, taking into account the diherent molar masses 
of Nd and La via the two-Debye function method^'^, 

( 2 ) 
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FIG. 2. (Color online) Log-log plot of the low-temperature dependence of the specific heat of the 
three magnetic states: the AFM (crosses), the FMrem (squares), and the FMfc (circles). 


with 


Therefore we get 




0.98. 


C{T) - = Cei(T) + a,p(T) + C^^,{T). 


( 3 ) 

( 4 ) 


To determine the Cm»A T) contribution for each magnetic phase, we can attempt to model 
the experimental data taking into account the different terms: C'ei(T) = 7T from free charge 
carriers, C'hyp(T) = AT~‘^ from the high-temperature limit of the Schottky anomaly due to 
the hyperfine splitting of the nuclear levels of the Nd^"*" ions, and C* m ag(T) from spin waves. 
The approach to study the last term was to consider the more general dispersion relations 
for the long-wavelength spin interactions. We examined the cases of AFM, FM, and type-A 
AFM (ferromagnetic layers antiferromagnetically coupled) states. The last is one of the 
proposed magnetic structures to occur below Tjv at zero applied magnetic field^. Because 
of the large magnetic anisotropy and the magnetoelastic effects present in this systen>^'^, we 
also took into account the possibility of the presence of a gap in the dispersion relation, A. 
In the low-temperature limit, the specific heat from each dispersion relation is found to be 


C'mag(r) = Be T 


T(12T2 + 6TA + A^) 

(5a) 

T"^(15T2 + 12TA + 4A2) 

(5b) 

(6T2 + 4TA +A^). 

(6c) 
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FIG. 3. (Color online) Zero-field low-temperature specific heat for both the AFM (diamonds) and 
FM (triangles) states, plotted in a log-log scale to remark the differences between both states. The 
lattice contribution to the specific heat has been subtracted to the experimental data. The lines 
are the best fits of Eq. (jj]) to the data. 


TABLE I. Results of the fitting of Eq. ([4]) to the specific heat data, for both the AFM and FM 
states. The units are mJ/(mol K™'"*“^) where m is the power of T corresponding to each coefficient, 
m = 1 for 7 , m = —2 for A, m = 2 for B. The number in parentheses is the statistical uncertainty 
in the last digit from the least-squares fitting procedure. 


state 

7 

A 

B 

A(K) 

AFM 

115(2) 

39.3(4) 

22.2(1) 

4.34(8) 

FM 

75(1) 

40.6(3) 

23.1(1) 

4.75(5) 


Equations (l5all . (l5bil and (|5B correspond, respectively, to the low-temperature magnetic 
contribution to the specific heat of the AFM, the FM, and the type-A AFM states. The 
usual expressions (T^, and T^, respectively) are recovered when A is set to zero. We 
fitted the experimental data to the expressions with and without gap and found that, nei¬ 
ther in the AFM nor in the FM phase, the “pure” AFM/FM contributions gave physically 
reasonable values for the parameters. On the contrary, in both phases the gapped type-A 
AFM contribution [Eq. (JScl)] was found to fit precisely. Figure [3] shows the specific heat for 
the AFM and FM states together with the best fits of Eq. (jl]) to the data. Table [I] lists the 
coefficients of all contributions. 
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The mixture of interactions and dimensionalities resulting in a contribution to Cmag 
has been also considered in magnetic structures where FM droplets are found in an AFM 
phasei^. In our system, nevertheless, the magnetization measurements indicate a saturated 
FM state, without evidence of AFM interactions. Despite we do not have a clear explanation 
to the presence of this term in the FM phase, the values for the parameters obtained fitting 
other contributions [Eq. fl5aD and (I5bp ] do not have any physical meaning. The larger gap 
obtained in the FM state corresponds with the larger internal magnetic field in this phase, as 
it will be shown below. The values of 7 for both phases are reasonable within the electronic 
contributions of rare earth intermetallics RsGe^ and do not need extra considerations. The 
reduction of the value in the FM phase with respect to the AFM can be attributed to a 
decrease in the density of states at the Fermi level, that would probably favour one of the 
electronic spin projections. 

The hyperfine contribution can not be omitted to fit completely the measured specific 
heat. The Schottky anomaly consists in a peak originated from the (de)population of discrete 
energy levels. In this case, these correspond to the hyperfine split nuclear levels of the Nd^^ 
ions. The Schottky anomaly can be approximated to A/T^ in its high-temperature limit, 
where A is related to the internal hyperfine magnetic field by the expressions^ 

( 0 ) 

Here I is the nuclear spin, fij is the nuclear magnetic moment, hfhyp is the internal mag¬ 
netic field at the Nd site, and the factor 5 is the number of moles of Nd per mole of 
NdsGes. Only two isotopes of Nd have nuclear spin different from zero (/ = 7/2), ^^^Nd 
and ^^®Nd with the natural abundances of 12.18% and 8.29%, whose nuclear magnetic mo¬ 
ments are p/ = —1.208/rAr and fij = — 0.744pjv respectivelySS. The hyperfine field values 
obtained are /io77hyp(AFM) = 272(2) T and /ro77hyp(FM) = 276(1) T. The energy splitting 
{AE = found is ~ 2.5 /leV for both phases. We can compare this value with 

the splitting of other Nd compounds studied with neutron scattering. Figure 0] plots the 
hyperfine energy splitting versus the saturated ionic magnetic moment of Nd for several 
Nd-based compounds^, along with the data point obtained in this work. The value used for 
the magnetic moment of the Nd ion corresponds to the one observed in the saturated FM 
stated, /iNd ~ 2/iB. It is remarkable that we have obtained approximately the same splitting 
for both magnetic phases, Ai7(AFM) = 2.50(2) /reV and Ai7(FM) = 2.53(1) peV, which 
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FIG. 4. (Color online) Energy splitting of inelastic neutron scattering signals in several Nd-based 
compounds (circles) as a function of the corresponding saturated ionic magnetic moment of Nd 
at low temperatures (adapted from Chatterji et al.—). The data point obtained in this work for 
NdsGes (square) is also shown. 

means that the average local magnetic moment per Nd ion is roughly the same in the two 
phases and corresponds to the value in the FM saturated state. Therefore, we may assert 
that the magnetic-field-induced AFM—)-FM transition simply flips the magnetic moments 
and preserves the value of /iNd- 

We will now proceed to investigate the temperature dependence of the magnetic contri¬ 
bution to the specific heat. Subtracting to the total measured specific heat the analytical 
functions of the electronic contribution, the Schottky anomaly and the phononic contribu¬ 
tion from the corrected LasGes data we obtain Gmag(7"). Fig. [5] shows the zero-field-cooled 
(ZFC) and remanent (rem) curves, where the latter was acquired as the sample was heated 
in zero applied magnetic field after it had been driven to the FM saturated state. The peak 
at 50 K indicates the temperature at which the AFM ordering takes place (T^r). The inflec¬ 
tion point observed around 26 K is in contradiction with the previously reported absence of 
any anomaly around this temperature, which was related to the possibility of the system to 
be in a spin glass state^i^^. 

From Gmag(F) we can compute the magnetic entropy as 



where it is assumed that the magnetic entropies of AFM and FM materials at zero tern- 
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FIG. 5. (Color online) Temperature dependence of the magnetic contribution to the specific heat 
obtained following the ZFC (diamonds) and rem (triangles) processes. A maximum at 50 K and an 
inflection point around 26 K are highlighted with vertical dotted lines. The lines joining the data 
points are guides to the eye. 
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FIG. 6. (Color online) Temperature dependence of the magnetic entropy obtained for the ZFC (dia¬ 
monds) and rem (triangles) processes. The horizontal solid line represents the value of 5i?ln(2J-|-l) 
for J = 9/2. The lines joining the data points are guides to the eye. 

perature are zero. Fig. E] shows how S'^ag attains the value of i?ln(2J + 1) expected for a 
paramagnet^ as the temperature grows above T^- In our case the entropy tends to 5i?ln 10, 
corresponding to 5 Nd^"*" free ions with J = 9/2. The actual value at which the obtained 
entropy tends is moderately smaller because the zero-held splitting due to the anisotropy 
could play a signihcant role even in the paramagnetic state. 
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FIG. 7. (Color online) Magnetic field dependence of the specific heat at 1.2 K, starting with the 
system in the AFM state (C ~ 0.2 J mol~^ K“^) and ending with the system in the FM state 
{C ~ 0.14 J mol~^ K“^). Two independent runs are plotted. The lines joining the data points are 
guides to the eye. 


Finally, the dependence of the specihc heat with the applied magnetic held at hxed tem¬ 
perature was studied to explore the magnetic-held-induced AFM—)-FM transition . The 
system was prepared following a ZFC process from T ^ down to 1.2 K, where the rel¬ 
ative difference between the specihc heat of the two phases has a maximum, and then the 
specihc heat was measured varying the held, from 0 to 40 kOe and back to 0 Oe. Fig. [7] 
shows how a large, abrupt, and irreversible change in C{H), associated with the AFM-FM 
transition, takes place between 16 and 17 kOe as the magnetic held ncreases. Two inde¬ 
pendent runs are plotted in the hgure showing the reproducibility of this transition. From 
the magnetic held-temperature phase diagram obtained from magnetization measurements 
in Ref. y one expects this transition to happen at much higher helds {H > 35 kOe). One 
explanation for this reduction of the held at which the spontaneous transition occurs is to 
consider it of thermally assisted origin. The large diherence in the thermal bath properties 
between the two experimental setups (MPMS® for magnetic measurements, versus PPMS® 
for specihc heat measurements) can strongly ahect how thermally-assisted abrupt transi¬ 
tions develop^i^i. We also see in the hgure that the decreasing dependence of the specihc 
heat with the increasing magnetic held, above H ^ 10 kOe, is consistent with the behaviour 
of a Omag term with a gap in the dispersion relation of the spin waves proportional to the 
applied magnetic held, A cx if. Nevertheless, a change in the sign of the slope is observed 
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in both states around 10 kOe, for which we do not have an explanation. A more detailed 
study considering also magnetoelastic effects could give a better description of the exact 
behaviour of the specihc heat as a function of the applied magnetic held. 


IV. CONCLUSIONS 

In summary, we have performed measurements of the specihc heat in the two magnetic 
phases of the system NdsGea. From the low-temperature data we have modeled the dif¬ 
ferent contributions to the specihc heat. A magnetic contribution is found in both the 
ferromagnetic (FM) and antiferromagnetic (AFM) phases. This term can be understood 
as a mixture of FM and AFM interactions in diherent dimensionalities. In the case of the 
AFM phase this term can be attributed to a type-A AFM, while in the case of the FM 
phase can be interpreted as a remanence of AFM interactions. The large magnetocrystalline 
anisotropy of NdsGes is evidenced by a gapped spin-wave spectrum in both phases. 

The average magnetic moment at low temperature in the two magnetic phases has been 
obtained by means of the specihc heat contribution of the hyperhne splitting of the nuclear 
moment of the Nd^"*" ions. The value of this magnitude is approximately 2/iB in both 
phases, which corresponds to the saturation value of the FM state at low temperature. 
Hence, we state that the magnetic-held-induced transition between both states corresponds 
to an irreversible spin-hip transition of the Nd ions. 

Finally, from the magnetic held dependence we observe that the held at which the sponta¬ 
neous transition takes place is remarkably smaller than the expected value from the magnetic 
held-temperature phase diagram. This is most likely due to the ehect of being in an envi¬ 
ronment with a smaller thermal coupling (PPMS® vs MPMS®), leading to a spontaneous 
ignition of a thermally assisted transition at smaller helds, probably by means of a magnetic 
dehagration process. 
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